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Abstract. We report about recent results on Dirac wave packets in the treatment of neu- 
trino flavor oscillation where the initial localization of a spinor state implies an interfer- 
ence between positive and negative energy components of mass-eigenstate wave packets. 
A satisfactory description of fermionic particles requires the use of the Dirac equation 
as evolution equation for the mass-eigenstates. In this context, a new flavor conversion 
formula can be obtained when the effects of chiral oscillation are taken into account. 
Our study leads to the conclusion that the fermionic nature of the particles, where 
chiral oscillations and the interference between positive and negative frequency com- 
ponents of mass-eigenstate wave packets are implicitly assumed, modifies the standard 
oscillation probability. Nevertheless, for ultra-relativistic particles and sharply peaked 
momentum distributions, we can analytically demonstrate that these modifications in- 
troduce correction factors proportional to ml^/Po which are practically un-detectable 
by any experimental analysis. 
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I. Introduction 

The Dirac formalism is useful and essential in keeping clear many of the conceptual aspects of quantum 
oscillation phenomena that naturally arise in a relativistic spin one-half particle theory. The quantum 
oscillation phenomena has stimulated the analysis of several theoretical approaches [1,2] on the flavor 
conversion formula which, sometimes, deserve a special attention because of carrying valuable physical 
information. The applicability of the standard plane wave treatment of oscillations by resorting to 
intermediate [3,4] and external [5,6] wave packet frameworks has been extensively questioned in the last 
years [1,7]. Although the standard oscillation formula [8,9] could give the correct result when properly 
interpreted, the plane wave approach implies a perfectly well-known energy-momentum and an infinite 
uncertainty on the space-time localization of the oscillating particle which leads to the destruction of 
the oscillating character [10]. The intermediate wave packet approach eliminates the most controversial 
points rising up with the plane wave formalism. Wave packets describing propagating mass-eigenstates 
guarantees the existence of a coherence length [10], avoids the ambiguous approximations in the 
plane wave derivation of the phase difference [7] and, under particular conditions of minimal slippage, 
recovers the oscillation probability given by the standard plane wave treatment [11]. Otherwise, a 
common argument against the intermediate wave packet formalism is that oscillating neutrinos are 
neither prepared nor observed [1]. Some authors suggest the calculation of a transition probability 
between the observable particles involved in the production and detection process in the so-called 
external wave packet approach [1,5]: the oscillating particle, described as an internal line of a Feynman 
diagram by a relativistic mixed scalar propagator, propagates between the source and target (external) 
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particles represented by wave packets. Anyway, it can be demonstrated [1] that the overlap function 
of the incoming and outgoing wave packets in the external wave packet model is mathematically 
equivalent to the wave function of the propagating mass-eigenstate in the intermediate wave packet 
formalism. However, the overlap function takes into account not only the properties of the source, but 
also of the detector. This is unusual for a wave packet interpretation and not satisfying for causality [1]. 
This point was clarified by Giunti [5] who solves this problem by proposing an improved version of 
the intermediate wave packet model where the wave packet of the oscillating particle is explicitly 
computed with field-theoretical methods in terms of external wave packets. In order concentrate the 
discussion on the Dirac equation properties that we intend to report in this manuscript, in this 
preliminary investigation, we avoid the field theoretical methods in detriment to a clearer treatment 
with intermediate wave packets which commonly simplifies the understanding of physical aspects going 
with the oscillation phenomena [7, 12]. 

Our final aim is the investigation of how the inclusion of chiral oscillation effects can modify the 
flavor conversion probability formula which was previously obtained by using fermionic instead of 
scalar particles, i. c. in treating the time evolution of the spinorial mass-eigenstate wave packets, 
we shall take into account the chiral nature of charged weak currents and the time evolution of the 
chiral operator with Dirac wave packets. To do it, we shall use the Dirac equation as the evolution 
equation for the mass-eigenstates. Before introducing the Dirac formalism, in section II we briefly 
review the analytic calculations [11] with the intermediate wave packet model for scalar particles [10]. 
In particular, a gaussian wave packet is chosen to describe the localization of the initial flavor state 
and to obtain an analytical expression for the flavor conversion probability. In section III, we shall 
recapitulate the Dirac formalism [13, 14] and show that a superposition of both positive and negative 
frequency solutions of the Dirac equation is often a necessary condition to correctly describe the 
time evolution of mass-eigenstate wave packets. The small modifications obtained in the context of a 
wave packet treatment of oscillation phenomena are (briefly) compared with quantum field theoretical 
calculations [1, 15, 16]. In section IV, we notice that the use of Dirac equation solutions allows us to 
observe the additional effect of chiral oscillation already introduced by De Leo and Rotelli [17]. As 
a natural extension, we show how to couple chiral to flavor oscillations in the intermediate wave 
packet framework. Finally, we give, for strictly peaked momentum distributions and ultra-rclativistic 
particles, an analytic expression for the coupled flavor and chiral conversion probability. We draw our 
conclusions in Section V. 



II. Scalar Oscillating Particles 

The main aspects of oscillation phenomena can be understood by studying the two flavor problem. In 
addition, substantial mathematical simplifications result from the assumption that the space depen- 
dence of wave functions is one-dimensional (z-axis). Therefore, we shall use these simplifications to 
calculate the oscillation probabilities. In this context, the time evolution of flavor wave packets can 
be described by 



where u a and up are flavor-eigenstates and u l and f 2 are mass-eigenstates. The probability of finding 
a flavor state up at the instant t is equal to the integrated squared modulus of the vp coefficient 



&(z, t) = (j>!(z, t) cos 6 v 1 + 4> 2 (z, t) s'm9 u 2 

= [<f>!(z, t) cos 2 9 + 4> 2 {z, t) sin 2 0] v a + [<f>i(z, t) — <j> 2 {z, t)] cosQsmOvp 
= </> a (z,t; 0)u a + </>p(z,t; 0)i/p, 



(1) 




+■ 



2 



{ 1 - SFO(t) } , 



(2) 



where SFO(t) represents the mass-eigenstate interference term given by 




+ 



(3) 
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Let us consider mass-eigenstate wave packets given at time t = by 



2 

7ra 5 



cxp 



exp [ip! z] 



where s = 1, 2. The wave functions which describe their time evolution are 

4>i{z,t) = J ^ tp(p m - pO exp[-i Eip^m^t + ip z z], 

(p, - p,) 2 a 2 



(4) 



(5) 



where 



E{p z ,m l ) — (j> 2 + mf) 2 and (p(p z — pj = (27ra 2 ) 4 exp 



In order to obtain the oscillation probability, we can calculate the interference term Sfo(£) by solving 
the following integral 



dp z 



fip, ~ Pi)<f(p* - V2) exp [-i AE(p z ) t] 



cxp 



{aAp) 2 



(0) 



where we have changed the z-integration into a /^-integration and introduced the quantities 

Ap = Pi-p 2 , Po = ^(Pi+p 2 ) and AE(jjS) = E{p zl m 1 ) E(p zl in 2 ). 

The oscillation term is bounded by the exponential function of a Ap at any instant of time. Under 
this condition we could never observe a pure flavor-eigenstate. Besides, oscillations are considerably 
suppressed if aAp > 1. A necessary condition to observe oscillations is that aAp <C 1. This constraint 
can also be expressed by Sp 3> Ap where Sp is the momentum uncertainty of the particle. The overlap 
between the momentum distributions is indeed relevant only for Sp 3> Ap. Consequently, without loss 
of generality, we can assume 



Sfo(*) = Re I J ^ V 2 {p z ~ Pa) exp H A£( P J f] 



(7) 



In litterature, this equation is often obtained by assuming two mass eigenstate wave packets 
described by the "same" momentum distribution centered around the average momentum p = p a . 
This simplifying hypothesis also guarantees instantaneous creation of a pure flavor eigenstate v a at 
t = [7]. In fact, for 0i(z,O) = <f> 2 (z,0) we get from Eq. {IJ 



b a (z,0,e)= — exp 



exp [ip a z] and <j>p(z, 0, 9) = 0. 



(8) 



In order to obtain an expression for ^(z, t) by analytically solving the integral in Eq. JSJ we firstly 
rewrite the energy E(p z ,m 1 ) as 



E(p z ,m i ) = E 1 



p: -Po 
E? 



= E 1 [1 + o~x (a, + 2v, 



(9) 



where 



7^ I 2 , 1\ i PO j P: - PO 

Ei = (m, +p ) 2 , v, = — and cr, = — - — 



The use of free gaussian wave packets is justified in non-relativistic quantum mechanics because the 
calculations can be carried out exactly for these particular functions. The reason lies in the fact 
that the frequency components of the mass-eigenstate wave packets, E(j> z ,m l ) = p 2 ,/2m 1 , modify 



4 



Alex E. Bernardini, Stefano De Leo: Flavor and Chiral Oscillations with Dirac Spinors 



the momentum distribution into "generalized" gaussian, easily integrated by well known methods of 
analysis. The term pi in E(p : ,,m 1 ) is then responsible for the variation in time of the width of the 
mass-cigcnstatc wave packets, the so-called spreading phenomenon. In relativistic quantum mechanics 
the frequency components of the mass-eigenstate wave packets, E(p z , mj = y/p 2 + m\ , do not permit 
an immediate analytic integration. This difficulty, however, may be remedied by assuming a sharply 
peaked momentum distribution, i. e. (ai^) -1 ~ a. <C 1. Meanwhile, the integral in Eq. (JjJJ can be 
analytically solved only if we consider terms up to order of in the series expansion. In this case, we 
can conveniently truncate the power series 



E(p zl m i ) = E t 



1 



(10) 



and get an analytic expression for the oscillation probability. The zero-order term in the previous 
expansion, E u gives the standard plane wave oscillation phase. The first-order term, p (Xi, will be 
responsible for the slippage due to the different group velocities of the mass-eigenstate wave packets 
and represents a linear correction to the standard oscillation phase [7]. Finally, the second-order 

2 

term, Tr^-of, which is a (quadratic) secondary correction will give the well-known spreading effects 
in the time propagation of the wave packet and will be also responsible for a new additional phase 
to be computed in the final calculation. In the case of gaussian momentum distributions for the 
mass-eigenstate wave packets, these terms can all be analytically quantified [11]. By substituting {TO} 
in Eq. (JSJ and changing the p 2 -integration into a cr 1 -intcgration, we obtain the explicit form of the 
mass-eigenstate wave packet time evolution, 



^(z, t) ?s (27ra 2 )-i cxp \—i{E 1 t — p z)] 

~2^ 



exp [—i (E l t — pa z)] exp 





a 2 E?af 






exp 


4 


exp 


—i 



ml t 
2iT 



■nal{t) 



<*?(*) 



i0.it, z) 



(11) 



where 



a,(i) = a 1 + 



Am* 
a 4 £? 



and Oiit, z) 



arctan 



2mH 



a 2 E? 



2ml t (z-Vity 



a 2 E? 



al(t) 



The time-dependent quantities a^t) and 9i(t,z) contain all the physically significant information [11] 
which arise from the second order term in the power series expansion 11U|) . By solving the integral J7J 
with the approximation @ and performing some mathematical manipulations, we obtain 



SFO(t) = BND(t) x Osc(i), 
where we have factored the time- vanishing bound of the interference term given by 

(Avty 



BND(i) = [1 + Sp 2 (t)] _1 exp 



2a 2 [1 + Sp 2 (£)]_ 

and the time-oscillating character of the flavor conversion formula given by 

Osc(t) = Re {exp [-iAEt - i©(f)]} = cos [AEt + Q(t)], 

where 



Sp(t) = - 2 A 
a 2 



with 



Av* 




P a 2 




and £ 




Po 










1 - 


3 + 









e(t) = 



SF 3 (t) 



- arctan fSp(t)l - ^t^t 
2 1 WJ 2p 2 [1 + SP 2 (t)] 



K 

E 2 



and 



E=y/E 1 E 2 . 



(12) 



(13) 



(14) 



(15) 



(16) 
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The time-dependent quantities Sp(£) and Q(t) carry the second order corrections and, consequently, 
the spreading effect to the oscillation probability formula. If AE <C E, the parameter p is limited by 

2 

the interval [1, —2] and it assumes the zero value when J^- « i. Therefore, by considering increasing 
values of p , from non-relativistic (NR) to ultra-relativistic (UR) propagation regimes, and fixing 

a f^ 2 , the time derivatives of Sp(i) and 6(i) have their signals inverted when -j^ reaches the value 
i. The slippage between the mass-eigenstate wave packets is quantified by the vanishing behavior 
of BND(i). In order to compare BND(t) with the correspondent function without the second order 
corrections (without spreading), 



BND ws (t) = cxp 



{Avty 

2a 2 



(17) 



we substitute Sp(t) given by the expression H14f) in Eq. and we obtain the ratio 



BND(t) 
BND ws (t) 



' AE t ' 
.o 2 E 2 



exp 



p 2 pi (AE t) 



(18) 



The NR limit is obtained by setting p 2 = 1 and p = in Eq. i|17|) . hi the same way, the UR 
limit is obtained by setting p 2 = 4 and p = E. In fact, the minimal influence due to second order 

2 

corrections occurs when J|r ~ ^ (p ~ 0). Returning to the exponential term of Eq. i|13|) . we observe 
that the oscillation amplitude is more relevant when Avt <C a. It characterizes the minimal slippage 
between the mass-eigenstate wave packets which occur when the complete spatial intersection between 
themselves starts to diminish during the time evolution. Anyway, under minimal slippage conditions, 

, , BND(t) -, 

we always have , , « 1. 

J BNDws(t) 

The oscillating function Osc(i) of the interference term SFO(t) differs from the standard oscillating 
term, cos[AEt], by the presence of the additional phase Q(t) which is essentially a second order 
correction. The modifications introduced by the additional phase 9(t) are discussed in Fig. [11] 
where we have compared the time-behavior of Osc(£) to cos [AE't] for different propagation regimes. 
The bound effective value assumed by 0(i) is determined by the vanishing behavior of BND(t). 

To illustrate this scalar flavor oscillation behavior, we plot both the curves representing Bnd(£) 
and 9(t) in Fig.[21[ll]. We note the phase slowly changing in the NR regime. The modulus of the phase 

2 

|0(i)| rapidly reaches its upper limit when > \ and, after a certain time, it continues to evolve 
approximately linearly in time. But, effectively, the oscillations rapidly vanishes. By superposing the 
effects of Bnd(£) in Fig. Eland the oscillating character Osc(t) expressed in Fig. ^ we immediately 
obtain the flavor oscillation probability which is explicitly given by 



I-II + Sp 2 ^)]" 3 exp 



(Avt) 2 



2a 2 [1 + Sp 2 (t)] 



cos[AE t + Q{t)} L(19) 



Obviously, the larger is the value of aE, the smaller are the wave packet effects. If it was suffi- 
ciently larger to not consider the second order corrections expressed in Eq. © , we could compute the 
oscillation probability with the leading corrections due to the slippage effect, 



sin 2 [26] 



1 — exp 



(Avt) 2 
2 a 2 



cos [AEt] 



(20) 



which corresponds to the same result obtained by [7]. By assuming an UR propagation regime with 
t f« L and Ei ~ p , under minimal slippage conditions (AvL <ti a), the above expression reproduces 
the standard plane wave result, 



Scalar ("a * ^jOi 



sin 2 [29] 
2 

sin 2 [29] 



1 - 



1 - 



(AV L) 
2 a 2 



1 — COS 



4po 



L 



sin 2 [26»] sin 2 



Am t 



COS 



[AE L] | 



(21) 



since we have assumed a E 3> 1 . 
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III. Dirac Formalism 

The results in the previous section have been obtained by considering scalar mass-eigenstates. Neu- 
trinos are, however, fermions. The time evolution of a spin one-half particle must be described by the 
Dirac equation. To introduce the fermionic character in the study of quantum oscillation phenomena, 
we shall use the Dirac equation as the evolution equation for the mass-eigenstates. The Eq. now 
becomes 

&(z, t) = ifj 1 (z 7 1) cosOv-l + ijj 2 (z, t) sin9v> 2 

= [ipi(z, t) cos 2 9 + ip 2 (z, t) sin 2 9] u a + [ipi(z, t) — ip 2 (z, t)] cos 9 sin 9 vp 

= ip a {z,t;6)v a +<ipi 3 (z,t;9)v l 3, (22) 

where ipi(z,t) satisfies the Dirac equation for a mass m^. The natural extension of Eq. (|SJ) reads 

ip a (z,0,e)=<l> a (z,0,6)w (23) 

where w is a constant spinor which satisfies the normalization condition w^w = 1. 



1 1 1.1. Dirac wave packets and the oscillation formula 

To describe the time evolution of mass-eigenstate Dirac wave packets, we could be inclined to superpose 
only positive frequency solutions of the Dirac equation. It seems, at first glance, a reasonable choice. 
However, when the initial state has the form given in Eq. (|23(1 . it is necessary to superpose both 
positive and negative frequency solutions of Dirac equation. In order to clear up this point, let us 
express the flavor state ip a {z, t, 9) in terms of 

f + °° dp 

ipi(z,t) = / — — exp [ip z z] {b s (p z , m t ) u s (p z , m t ) exp [-iE(p z ,m l )t] 

-» 77 i=l,2 

+ ^'(-p^m,) v s {-p z , mO exp [+iE(p„m 1 )t]}. (24) 

At time t = the mass-eigenstate wave functions satisfy i/j^ZjO) = tp 2 (z,0) (this guarantees the 
instantaneous creation of a pure flavor-eigcnstate u a as we have appointed in section II) . The Fourier 
transform of il>i{z 1 0) is 

[b s {p z ,m 1 )u s {p z ,m 1 ) + d^i-p^^m^v^-p^m,)]. (25) 

i=l,2 

By observing that the Fourier transform of (f> a {z, 0, 9) is given by <p(p z — p ) (see Eq. (|SJ|), we imme- 
diately obtain the Fourier transform of ipa{z, 0, 9), 

(p(p z -p )w= ^2 [b\p z ,m l )u\p zl m l ) + d s *(-p z ,m 1 )v s (-p z ,m 1 )]. (26) 

i=l,2 

Using the orthogonality properties of Dirac spinors, we find [14] 

6 s (p z ,TO!) = tp(p, -p )u st (p,,m 1 )w, 
^""(-p^m,) = tp(p z -Pojv^i-p^mjw. (27) 

These coefficients carry an important physical information. For any initial state which has the form 
given in Eq. I|23() . the negative frequency solution coefficient d s * (-p z ,m 1 ) necessarily provides a non- 
null contribution to the time evolving wave packet. This obliges us to take the complete set of Dirac 
equation solutions to construct the wave packet. Only if we consider a momentum distribution given 
by a delta function (plane wave limit) and suppose an initial spinor w being a positive energy mass- 
eigenstate with momentum p 0l the contribution due to the negative frequency solutions d s * (~p z ,m 1 ) 
will be null. 
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Having introduced the Dirac wave packet prescription, we are now in a position to calculate the 
flavor conversion formula. The following calculations do not depend on the gamma matrix represen- 
tation. By substituting the coefficients given by Eq. I|27|) in Eq. H24J1 and using the well-known spinor 
properties [14], 



we obtain 



/dp \ «7° (j 3 p~ + m-i) I 
— - ip(p z - pa) exp [ip z z] < cos [Eip^m^t] )- — - — - — - sin [E(p z: m,)i] } w.(29) 
oo 2lr { E{p z ,m l ) J 

By simple mathematical manipulations, the new interference oscillating term will be written as 
f + °° dp 

DFO(i) = / ip 2 (p z -po) { [1 - f(p z ,m u m 2 )} cos[e_{p z ,m 1 ,m a )t] + 

f(p z ,m u m 2 ) cos[e + (p z ,m 1 ,m 2 ) t}} (30) 

where 

ft \ E(p zl m 1 )E(p z ,m 2 ) -p 2 z - mi m 2 

f[p„m 1 ,m a ) = —— — r and e±{p z ,m 1 ,m 2 ) = E(p z ,m 1 ) ± E{p z ,m 2 ) . 

2E(p z ,m 1 )E(p z ,m 2 ) 

The time-independent term f(p zl m u m 2 ) deserves some comments. It has a minimum at p z = and 
two maxima at p z = ±^/mim 2 . It goes rapidly to zero when p z 3> m\.2 (ultra-relativistic limit) as well 
as when p„ -C /fii,2 (non-relativistic limit). It means that when we consider a momentum distribution 
sharply peaked around p„ ^ mi, 2 or p <ti mi,2 the corrections introduced by f{p z ,m l ,m 2 ) are 
negligible. The maximum value of f(p z ,m ll m 2 ) is 



f ma x{p,,m ll m 2 ) = \ ^ — - (31) 

z m 1 + m 2 

which vanishes in the limit m 1 = m 2 . The effects introduced by f(p z ,m 1 ,m 2 ) are relevant only when 
Am ?s m 1 3> m 2 . Meanwhile, what is interesting about the result in Eq. is that it was obtained 
without any assumption on the initial spinor w. Otherwise, the initial spinor carries some fundamen- 
tal physical information about the created state. And this could be relevant in the study of chiral 
oscillations [17] where the initial state plays a fundamental role. In comparison with the standard 
treatment of neutrino oscillations done by using scalar wave packets, where the interference term 
SFO(i) is given by Eq. J7J) with AE(p z ) = e_(j) z ,m l ,m 2 ), we note in Dfo(£) two additional terms. In 
the first one, the standard oscillating term cos [e_ (p z , m, , m 2 ) t] , which arises from the interference be- 
tween mass-eigenstate components of equal sign frequencies, is multiplied by a new factor obtained by 
the products (p z ,m 1 ) u(p z ,m 2 ), (-p z ,m 1 ) v(-p z ,m 2 ) and h.c. The second one is a new oscillating 
term, cos [e + (p z7 m t , m 2 ) t], which comes from the interference between mass-eigenstate components 
of positive and negative frequencies. The factor multiplying such an additional oscillating term is 
obtained by the products u^(p z ,m 1 )v(-p z ,m 2 ), (-p z: m 1 ) u(p z ,m 2 ) and h.c. The new oscillations 
have very high frequencies. Such a peculiar oscillating behavior is similar to the phenomenon referred 
to as Zitterbewegung. In atomic physics, the electron exhibits this violent quantum fluctuation in the 
position and becomes sensitive to an effective potential which explains the Darwin term in the hydro- 
gen atom [18]. We shall sec later that, at the instant of creation, such rapid oscillations introduce a 
small modification in the oscillation formula. 

Returning to the starting point, if we had postulated a wave packet made up exclusively of pos- 
itive frequency plane-wave solutions, the oscillation term cos [e + (p z , m l7 m 2 ) t] would have vanished. 
It reinforces the argument that, in constructing Dirac wave packets, we cannot simply forget the 
contributions of negative frequency components. 
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1 1 1.2. First order modifications to the oscillation formula 

A more satisfactory interpretation of the modifications introduced by the Dirac formalism is given 
when we explicitly calculate DFO(i). By considering the energy E(p z7 m l ) expansion up to the second 
order terms in Eq. I|10|) . we include an analysis of spreading effects. In this preliminary study, we are, 
however, interested only to first order corrections. Thus, we approximate the frequency components 

by 

E(p z , mi ) » E, + v, (p, -po) . (32) 
As a consequence of this approximation, we get 



f(p z ,m u m 2 ) 



1 - v lV2 1 + 



m 1 m 2 
Pi 



ViV 2 



+ <) 1 + 



P* - Po 
Po 



and 



e±(p J8 ,m 1 ,m 2 ) w E t ±E 2 + (v x ±v 2 ) (p z -p ) 



(33) 



(34) 



For UR particles (m, <C p ), we can also use the following expression for the central energy values 
(Ei) and the group velocities (v^ of the mass-eigenstate wave packets, 



E 1 w p„ 



2p„ 



and v,« 1 



This implies 



/( R , mi , m2 ),(^)7l-2^ 
\2p J V Po 



e + (p z ,m 1 ,m 2 ) « 2p 
Am 



1 



m\ + m 2 2 p z — p 



e_(p z ,m 1 ,m 2 ) 



2 Po 



1 



4p 2 
V* - Po 

Po 



Po 



1 - 



ml + ml 
4p 2 



where (Am) 2 = (m x — m 2 ) 2 is different from Am 2 = ml — ml which appears in the standard oscillation 
phase. Finally, by simple algebraic manipulations and after gaussian integrations, we find 



Dfo(<) « exp 



/ Am 2 t V 
V 2\/2ap2 J 



1 



Am 
2po 



Am 2 
2po 



Am 

2pu 



Ar, 



'Pa 



-t sin 



Am 2 
2po 



exp 



2a 2 



( 2 -^Y (^) 2 {cos[ Po t(2 + ^) 



+ 1» ( 2 - ,J ^r) sin [** ( 2 + ^)] }■ ( 35 ) 



As we have already noticed, the oscillating functions going with the second exponential function in 
Eq. (|35fl arise from the interference between positive and negative frequency solutions of the Dirac 
equation. It produces very high frequency oscillations which is similar to the quoted phenomenon of 
Zitterbewegung [18]. The oscillation length which characterizes the very high frequency oscillations is 
given by L^ S ^! F m Obviously, L^^ F is much smaller than the standard oscillation length given 

by Lq*c = . It means that the propagating particle exhibits a violent quantum fluctuation of its 
flavor quantum number around a flavor average value which oscillates with Lq'c • Meanwhile, except 
at times t ~ 0, it provides a practically null contribution to the oscillation probability. To explain 
such a statement, let us suppose that an experimental measurement takes place after a time t ~ L 
for UR particles. The observability conditions impose that the propagation distance L must be larger 
than the wave packet localization a. Since the (second) exponential function vanishes when L»a, for 
measurable distances, the effective flavor conversion formula will not contain such very high frequency 
oscillation terms, and can be written as 



•[26] 



1 — exp 



' Am 2 L \ ' 



1 



( A™l\ ' 
\ 2po j 



2po 
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For distances which are restrict to the interval a <C L <C a 2 ^^ we observe the minimal slippage 
between the wave packets. In this case, we could suddenly approximate the oscillation probability to 



sin 2 [20] 



1 - 



I _ ( Am 2 L V 



1 - 



2po 



COS 



2po 



(37) 



however, we reemphasize that it is not valid for T w L ~ when the rapid oscillations are still relevant 
(L < a). By comparing the result of Eq. I|37|l with the scalar oscillation probability of Eq. i|12|) . we 

notice a deviation of the order (^^fj that appears as an additional coefficient of the cosine function. 
It is not relevant in the UR limit as we have noticed after studying the function f(p z ,m 1 , m 2 ). 



1 1 1.3. A brief extension to quantum field treatment 

Now we try to establish a tenuous correspondence between our results and the quantum field theory 
(QFT) treatment. It was extensively demonstrated in the literature [5,6,16] that the oscillating particle 
cannot be treated in isolation. The oscillation process must be considered globally: the oscillating states 
become intermediate states, not directly observed, which propagate between a source and a detector. 
This idea can be implemented in QFT when the intermediate oscillating states are represented by 
internal lines of Feynman diagrams and the interacting particles at source/detector are described 
by external wave packets [1,16]. In this context, let us consider the weak flavor-changing processes 
occurring through the intermediate propagation of a neutrino, 

Pi — > Pf + ct + v a (oscillation) vp + Di — > (3 + Dp (38) 

where pi and pp (Dj and Dp) are respectively the initial and final production (detection) particles. 
The amplitude for the process is represented by 

A = (p F ,D F |T(cxp [-i J cte 4 W/]) - l| pi,Dr) (39) 

where TLi is the interaction Hamiltonian for the intermediate particle and T is the time ordering op- 
erator. After some mathematical manipulations [1], this amplitude can be represented by the integral 

/dE d 3 T) 
-^f F(E, p) G(E, p, t D , t P ) cxp [i p ■ (xj, - x P )] (40) 

where the function F(E,p) represents the overlap of the incoming and outgoing wave packets, both 
at the source and at the detector, and the Green function in the momentum space, G(E,p,tp,tp), 
represents the fermion propagator which carries the information of the oscillation process. The overlap 
function is independent of production and detection times and positions (tp, tp,, xp, X£>) and depends 
on the the directions of incoming and outgoing momenta. In certain way, the physical conditions of 
source and detector, in terms of time and space intervals, are better defined in this framework than 
in the intermediate wave packet framework. Anyway, to understand the oscillation process we must 
turn back to the definition of mixing in quantum mechanics. It is similar in field theory, except that 
it applies to fields, not to physical states. This difference allows to bypass the problems arising in the 
definition of flavor and mass bases [1]. In one-dimensional spatial coordinates, the mixing is illustrated 
by the unitary transformation 

^(z, t; 9) = £-'(0; t) ^(z, t) G{9- 1) (41) 

as the result of the noncoincidence of the flavor basis (a = a, (3) and the mass basis (i = 1, 2). The 
Eq. H41fl gives the Eq. when the generator of mixing transformations Q(8; t) is given by 

gtfi; t) = cxp[6> J dz ^(z, t)i) 2 {z, t) - tp 2 {z, t)Vi(«, t)]. (42) 

By taking the one-dimensional representation of Eq. H40(l . the propagator G(E,p z ,tp, tp) can also be 
written in the flavor basis as 



G a0 (9; E,p.,T) = g-'ie- 1) G(E, Pz7 T) G(6; t) = g~\0; t) G(E, Pz ,t D ,t P ) G(0; t) (43) 
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with T = to — tp. 

In particular, by following the Blasone and Vitiello (BV) prescription [15,19], the definition of a 
Fock space of weak cigenstates becomes possible and a nonperturbative flavor oscillation amplitude 
can be derived. In this case, the complete Lagrangian (density) is split in a propagation Lagrangian, 

C p (z, t) = ijj^z, t) (ifi - m x ) ip ± (z, t) + Tp 2 (z, t) (ifi - m 2 ) %p 2 (z, t), (44) 

and an interaction Lagrangian 

d(z, t) = ip a (z, t] 0) {ifi - m a ) ^ Q (z, t] 9) + ipp(z, t; 0) {ifi - mp) ippiz, t; 9) 

- m a p ($ a {z, t; 9)ip p (z, t; 9) + $p(z, t; 9)^ a (z, t; 9)) , (45) 

where 

m a(/3) = m i(2) cos2 $ + m 2(i) sin 2 and m a p = (m, — m 2 ) cos 9 sin 9. 

In general, the two subsets of the Lagrangian can be distinguished if there is a flavor transformation 
which is a symmetry of Ci{z, t) but not of C p (z, t). Particle mixing occurs if the propagator built from 
C p (z,t), and representing the creation of a particle of flavor a at point z and the annihilation of a 
particle of flavor (3 at point z', is not diagonal, i.e. not zero for (3 = a. The free fields ipi(z,t) can 
be quantized in the usual way by rewriting the momentum distributions b s (p z ^mi) and d s *(-p z ,m 1 ) 
in Eq. (|24J) as creation and annihilation operators B s (p z ,m 1 ) and D s ^{^p z , mj. The interacting fields 
are then given by 

VvCM) = j + Z ff ex PM E W(p.it)<(p.;t)+^(-p.;i)«;(-p.;i)} (46) 

i=l,2 

where the new flavor creation and annihilation operators which satisfy canonical anticommutation 
relations are defined by means of Bogoliubov transformations [19] as 

SJ(p,;t) =g- 1 (9;t)B s ( P :,m 1 )g(9;t) and D s a (- Pz ; t) = Q-\0; t) D s {-p z , m.) G(9; t). 

By following the BV prescription [15], which takes into account the above definitions, it was demon- 
strated [20] that the flavor conversion formula can be written as 

P(u a -> vp-t) = \{B s p(p ;t), B s a (p ;t)}\ 2 + \{D s (-p o ;t), B s a (p ;t)}\ 2 (47) 

which is calculated without considering the localization conditions imposed by wave packets, i. e. by 
assuming that p z ~ p . When the explicit form of the flavor annihilation and creation operators are 
substituted in Eq. (07J), it was also demonstrated [19] that the flavor oscillation formula becomes 

P(v a -> up;t) = sm 2 [26] {[1 - f(p ,m u m 2 )} cos[e_(p , m u m 2 ) t] 

+ f(p ,m 1 ,m 2 ) cos[e + (p ,m u m 2 )t}} 

- sin 2 [29] { [l - (^) 2 ] sin 2 [*£t] + sin 2 [ Po t (l + ^)] } (48) 

where the last approximation takes place in the relativistic limit p 3> ^/m 1 m 2 . After some simple 
mathematical manipulations, the Eq. (|48|l gives exactly the oscillation probability P Dirac (i^ a — > vp\V) 
calculated from Eq. I)35|l when it is assumed that the wave packet width a tends to infinity and t ~ L. 

This new oscillation formula tends to the standard one (ED in the UR limit. If the mass cigenstates 
were nearly degenerate, we could have focused on the case of a nonrelativistic oscillating particle hav- 
ing very distinct mass eigenstates. Under these conditions, the quantum theory of measurement says 
that interference vanishes. Therefore, as we have already appointed, the effects are, under realistic 
conditions, far from observable. Besides, in spite of working on a QFT framework, the lack of observ- 
ability conditions must be overcome by implementing external wave packets, i. e. by calculating the 
explicit form of Eq. I|40|l for fermions. Such a procedure was applied by Beuthc for scalar particles [1] 
and, in a very particular analysis, with basis on the BV calculations and on our intermediate wave 
packet results, it could be extended to the fermionic case. 
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IV. Flavor and Chiral Oscillations 



In treating the time evolution of the spinorial mass-eigenstate wave packets in the previous section, 
we have overlooked an important feature. We have completely disregarded the chiral nature of charged 
weak currents and the time evolution of the chiral operator. In the following, we aim to investigate if 
(and eventually how) the flavor oscillation formula could be modified by this additional effect. 

It is well known that from the Hciscnberg equation, we can immediately determine whether or 
not a given observable is a constant of the motion. If neutrinos have mass, the operator 7 s does not 
commute with the mass-eigenstate Hamiltonians. This means that for massive neutrinos chirality is not 
a constant of the motion. Observing that neutrinos with positive chirality are decoupled from charged 
weak currents, this additional effect cannot be ignored. We have already seen that localized states 
contains, in general, plane-wave components of negative and positive frequencies. As an immediate 
consequence of this, the interference between positive and negative frequencies, responsible for the 
additional oscillatory term in DFO(i), will also imply an oscillation in the average of chirality. Thus, 
the use of Dirac equation as evolution equation for neutrino mass-eigenstate wave packets leads to 
an oscillation formula containing both "flavor-appearance" (neutrinos of a flavor not present in the 
original source) and "chiral-disappearance" (neutrinos of wrong chirality) probabilities. 

We obtain the Dirac flavor and chiral oscillation probability formula in the same way as we have 
obtained the Eq. (|3fc>|> . By assuming that the normalizable mass-eigenstate wave functions ip^ 2 (z,t) 
are created at time t = as a — 1 chiral eigenstate, we can write 




dzipl(z,t)'f'ip i (z,t) I = 
1 - f(p t ,rrh,rrh) - 



dp z 
~2^~ 



V 2 {Pz -Po) 



E(p z ,m 1 )E(p^,m 1 ) 
m,m, 

I - — 

E(p z ,m 1 )E(p z ,m l ) 



COS [e_ (j>. , TOi , TOi ) t] 

cos{e + (p z ,m l ,m 1 )t] )■ = 1,2) .(49) 



From this integral, it is readily seen that an initial —1 chiral mass-eigenstate will evolve with time 
changing its chirality. Once we know the time evolution of the chiral operator, we are able to construct 
an effective oscillation probability which takes into account both flavor and chiral conversion effects, 
i.e. 



^,1^^^)= / dz\^ L (z,t;6)\ 




dz\iP ih {z,t)\< 



[DCO(t) - DFCO(f)] 



Re 



dzif)[ L (z,t) ^ 2 ^(z,t) 



(50) 



As done in the previous section, the terms Dco(i) and DFCO(t) can be rewritten by using a p z 
integration, 



1 2 [ + °° dp 

Dco(t) = - ^2 / g - ^ ^ 2 (P- ~Po) { 1 ~ c^m^mO + c^^m^TO,) cos[2 E(p z , m t ) t}} 



m; 



in; 



2ir ^(p, P'^ \ 2E 2 ( Pz , mi ) 2E 2 (p z , mi ) 



cos[2 Eip^mjt] \ (51) 



and 



DFCO(t) 



7T~ ^(P* ~Po) { [1 _ c(p z ,m 1 ,m 2 )} cos[e _(p z ,m 1 ,m 2 )t] 

Z 7T 



+ c{p z ,m 1 ,m 2 ) cos[e + (p z ,m u m 2 )t]} , (52) 
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where 



zip^m^rrij) = f{p z ,m 1 ,m J ) + 



2E{p % ,m 1 )E{p :s ,m j ) 



The functions c{p z ,m 1 ,m :i ) have a common maximum at p z = which, contrary to what happened 
for f(p e ,m 1 ,m 2 ), do not depend on the mass values, 

c max (0, m u m 3 ) = \ , 

and, following the same asymptotic behavior of f(p z , m l , m 2 ), go rapidly to zero for p z 3> m lt2 - As a 
consequence of the first order approximation (1.''>2|I , we get 



1 - VjVj + v,Vj (v 2 + v 2 - 2) 



Pz-P0 

Po 



which gives 



m 2 + to 2 
4p 2 



1 - 2 



P z - Po 

Po 



in the UR approximation. By substituting c(p z ,m l ,m j ) in the above integrations (j51l52() and after 
some algebraic manipulations, we explicitly calculate the terms Dco(t) and DFCO(t), 



Dco(t) « 1 - fife + exp 
X < cos 



Po 



2 p 2 -m\ 



*^-% m1 1 sin 

a 2 Po 



2 Pq +m 
Po 



i /'o 



exp 



2 Po -™S ^ 



i Po 



J cos \^pS±A t ] + 4p °-\ ml t 

I L Po J a Po 



sm 



2 p +m 2 
po 



! ']}} 
}}■ 



Dfco(£) « exp 



2\/2a 



Po 







cos 






L 1 4p 2 




_ 2p 



— a — i — ~o — r £ sm 

4 Po a Po 



2po 



- exp 



/ 4p 2 -m 2 -r; 
V 2V2ap 2 



„ 1 rj H — — — j-T — - i sm 

2p a 2 pg 



4p (] +m 1 +m 2 
2po 



(53) 



(54) 



Again, in the hypothesis of minimal slippage between the mass-eigenstate wave packets (AvL <C a), 
and for long distance between source and detector [L 3> a), i.e. 

L 



1 <C - < -f^—, 
a A m 2 

the standard flavor oscillation probability is reproduced. In fact, 



P 



'[26] 



1 



m 1 +m 2 



sin 2 [261] f 
~ 2 \ 

= sin 2 [20] sin 2 



1 — cos 



{'-[ 



- 1- 



2\/2ap 2 



COS 



Am 1 
2po 



I 



Am r 



(55) 



V. Conclusions 

In order to quantify some subtle changes which appear in the standard flavor oscillation probability 
[9] due to chiral oscillations coupled to the flavor conversion mechanism of free propagating wave 
packets, we have reported about some recent results on the study of flavor oscillation with Dirac 
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wave packets [11]. By taking into account the spinorial form of neutrino wave functions and imposing 
an initial constraint where a pure fiavor-eigenstate is created at t = 0, for a constant spinor w, 
it is possible to calculate the contribution of positive and negative frequency solutions of the Dirac 
equation to the wave packet propagation and, finally, to obtain the oscillation probability. Particularly 
we have noticed a term of very high oscillation frequency depending on the sum of energies in the new 
oscillation probability formula. In addition, the spinorial form of the wave functions and their chiral 
oscillating character subtly modify the coefficients of the oscillating terms in this flavor conversion 
formula. To describe the time evolution of the mass-eigenstates, we have assumed an initial gaussian 
localization and performed integrations by considering a strictly peaked momentum distribution. 
Under the particular assumption of UR particles, we have been able to obtain an analytic expression 
for the coupled chiral and flavor conversion formula. In case of Dirac wave packets, these modifications 
introduce correction factors which arc negligible in the UR limit. We have confirmed that the fermionic 
character of the particles modify the standard oscillation probability which is previously obtained by 
implicitly assuming a scalar nature of the mass-eigenstates. 

However, we know the necessity of a more sophisticated approach is understood. It involves a 
field-theoretical treatment. Derivations of the oscillation formula resorting to field-theoretical meth- 
ods are not very popular. They are thought to be very complicated and the existing quantum field 
computations of the oscillation formula do not agree in all respects [1]. The Blasone and Vitiello 
(BV) model [2,15] to neutrino/particle mixing and oscillations seems to be the most distinguished 
trying to this aim. They have attempt to define a Fock space of weak eigenstates and to derive a 
nonperturbative oscillation formula. Flavor creation and annihilation operators, satisfying canonical 
(anti)comutation relations, are defined by means of Bogoliubov transformations. As a result, new 
oscillation formulas arc obtained for fcrmions and bosons, with the oscillation frequency depending 
not only on the difference but also on the sum of the energies of the different mass-eigenstates. Mean- 
while, the prescription of oscillating neutrinos as Dirac spinors was not yet completely and accurately 
described in a quantum field formalism. With Dirac wave packets, the flavor conversion formula can 
be reproduced [13] with the same mathematical structure as those obtained in the BV model [2, 15]. 
Moreover each new effect present in the oscillation formula can be separately quantified. 

In fact, the quantum-mechanical treatment which associates Dirac wave packets with the propagat- 
ing mass eigenstates is rich in physical insights which were extensively studied in this paper. Besides 
the review of analytical calculations done with gaussian wave packets for scalar [11] and fermionic [13] 
particles, the main conceptual aspect arising from the formalism with Dirac wave packets leads to 
the study of chiral oscillations. In the standard model flavor-changing interactions, neutrinos with 
positive chirality are decoupled from the neutrino absorbing charged weak currents [17]. A state with 
left-handed helicity can be approximated by a state with negative chirality in the UR limit. Once we 
have assumed the interactions at the source and detector are chiral only the component with negative 
chirality contributes to the propagation. In this case, we are obliged to consider chiral coupled to 
flavor oscillations in order to compute the modifications to the standard flavor conversion formula. In 
fact, when chiral oscillations are taken into account, these modifications introduce correction factors 
proportional to 2 /p 2 which are, however, practically un-detectable by the current experimental 
analysis. It leads to the conclusion that, in spite of often being criticized, the standard flavor oscilla- 
tion formula resorting to the plane wave derivation can be reconsidered when properly interpreted, but 
a satisfactory description of fermionic (spin one-half) particles requires the use of the Dirac equation 
as evolution equation for the mass-eigenstates. 
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Figures 
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Fig. 1. The time-behavior of Osc(t) compared with the standard plane wave oscillation given by cos [AEt] 
for different propagation regimes. The additional phase O(i) changes the oscillating character after some time 
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of propagation. The maximal deviation occurs for -|% ~ | • 
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Fig. 2. The time-behavior of the additional phase Q(t). The values assumed by 0(t) are effective while the 
interference term does not vanish. In the upper box we can observe the behavior of Bnd(£) which determines the 

2 

limit values effectively assumed by Q(t) for each propagation regime. For relativistic regimes with ^ > 5, the 
function 0(t) rapidly reaches its lower limit as we can observe in the small box above. We have used aE — 10. 



